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We study numerically the localization properties of two-dimensional electrons in a weak perpen- 
dicular magnetic field. For this purpose we construct weakly chiral network models on the square 
and triangular lattices. The prime idea is to separate in space the regions with phase action of 
magnetic field, where it affects interference in course of multiple disorder scattering, and the re- 
gions with orbital action of magnetic field, where it bends electron trajectories. In our models, the 
disorder mixes counter-propagating channels on the links, while scattering matrices at the nodes 
describe exclusively the bending of electron trajectories. By artificially introducing a strong spread 
in the scattering strengths on the links (but keeping the average strength constant), we eliminate 
the interference and reduce the electron propagation over a network to a classical percolation prob- 
lem. In this limit we establish the form of the disorder - magnetic field phase diagram. This 
diagram contains the regions with and without edge states, i.e. the regions with zero and quantized 
Hall conductivities. Taking into account that, for a given disorder, the scattering strength scales 
as inverse electron energy, we find agreement of our phase diagram with levitation scenario: energy 
separating the Anderson and quantum Hall insulating phases floats up to infinity upon decreasing 
magnetic fleld. From numerical study, based on the analysis of quantum transmission of the network 
with random phases on the links, we conclude that the positions of the weak-field quantum Hall 
transitions on the phase diagram are very close to our classical-percolation results. We checked 
that, in accord with the Pruisken theory, presence or absence of time reversal symmetry on the Unks 
has no effect on the line of delocalization transitions. We also find that floating up of delocalized 
states in energy is accompanied by doubling of the critical exponent of the localization radius. We 
establish the origin of this doubling within classical-percolation analysis. 

PACS numbers: 72.15.Rn; 73.20.Fz; 73.43.-f 



I. INTRODUCTION 



Levitation scenario 



Scaling theory of localization « predicts that evolution 
with size, L, of the conductivity, cr, (in the units of 
e^/27r?i) of a 2D sample is governed only by the value 
of (T, regardless of the type of disorder in the sample, i.e., 
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Together with initial condition, cr|^^; = ctq = kpl, where 



kf. is the Fermi momentum, and I is the transport mean 
free path, Eq. ([T]) suggests that, in zero magnetic field, 
where fi{a) = —2/{TTa), localization radius of electron 
states is given by 
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With increasing magnetic field, when /3(fT) crosses from 
the orthogonal to the unitary form, /3„(cr) = — l/(7rcr)^, 
the localization radius crosses over from to given 
by 



In 
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Crossover takes place when weak localization is sup- 
pressed, i.e., when the magnetic flux through the vector 
area spanned by the electron travelling diffusively over 




2 1/V. 



FIG. 1: (Color online) Left: Energy position of two lowest 
delocalized states, Eq and as a function of magnetic fleld, 
LUc, as predicted in Ref. ^. At a given E, delocalization tran- 
sition occurs at low fleld (point L) and at high fleld (point 
H). Right: Same dependencies, replotted in the axis inverse 
flUing factor vs. zero-field resistivity, constitute a part of 
the global phase diagram Ref. [^. Depending on p°^, the 
system undergoes a sequence of transitions, — >■ 1 — !■ or 
1 ^ 2 ^ 1 ^ 0. 



the area ^ is of the order of the flux quantum. With 
vector area being ~ ^o^, we get the following estimate for 
crossover magnetic field 
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where uj^ is the cyclotron frequency, and r is the scatter- 
ing time. 

The origin of the crossover Eq. (jH) is that the paths, 



which interfere in a zero field, acquire field-induced ran- 
dom Aharonov-Bohm phases. Justification for consider- 
ing exclusively the phase action of magnetic field is that 
in high-mobility samples with kpl 3> 1 the crossover field 
is so weak that its orbital action can be neglected. It is 
a very delicate fact that, after the crossover, this orbital 
action causes a drastic change of the eigenstates even 
for classically weak magnetic fields, lUcT <^ 1. This con- 
clusion was drawn by Khmelinitskii^ from the analysis 
of the renormalization group Rows^. Originally, these 
equations were derived to describe quantization of the 
Hall conductivity in a strong-field limit, uJcT 3> 1, 
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where Uxx and Uxy are, respectively, the diagonal and 
non-diagonal components of the conductivity tensor, and 
I? is a dimensionless constant. First term of Eq. ([5]) is 
the same as in Eq. ^ with unitary /3(cr). It originates 
from interference: two paths corresponding to the same 
scatterers but different sequences of scattering events in- 
terfere even in the presence of Aharonov-Bohm phases. 
Second term reflects the orbital action of magnetic field 
(Lorentz force); by curving electron trajectories it tends 
to destroy the interference. Quantum Hall transition be- 
tween cyxy\]^^^ = n and (^xy\^^^ = n A- \ takes place 
when the "phase" and "orbital" terms compensate each 
other. Khmelnitskii's treatment^, is equivalent to solving 
Eqs. (HI, (|6|) together with classical Drude initial condi- 
tion. 
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which yields the positions of delocalized states 
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As shown in Fig. [T] for n — 0, 1, the high-field part, 
ujcT ^ 1, of En follows the centers of Landau levels, 
while the low-field part, i?„T ^ {n + 1/2){u!ct)~^ , "floats 
up" as LOcT — > 0. Such a behavior of critical values of 
kpl — EpT in vanishing field is usually called levitation 
of delocalized states^. More specifically, it is expected 
that, upon increasing magnetic field above the crossover 
Eq. (|4]), localization radius, ^(wc), diverges in the vicinity 
of discrete values WcT = {n + l/2)(fcFZ)~^, changing from 
unitary to infinity and returning back to S^u- Recast- 
ing Eq. ([S]) into the dependence, pjjj, — (Tq^^, versus the 
inverse filhng factor, = uJc/{2Ep), yields a system of 
semicircles 
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shown in Fig. [T] inset, which is a part of the global phase 
diagram^. This diagram suggests that for high enough 




FIG. 2: In non-chiral network model, confining potential re- 
stricts electron motion to the links. All the scattering ma- 
trices S at the nodes are the same. Positional disorder is 
emulated by randomness of phases acquired on the links. De- 
pending on symmetry class, the phases between the points 1 
and 2 are either the same or difi'erent. 



i.e., for strong disorder, the resistance pxxi^c) grows 
monotonically. Upon decreasing below a certain 
threshold value. Fig. [I](e.g., by applying the gate volt- 
age) Pxxi^c) exhibits two quantum Hall peaks. For p'^^ 
smaller than the second threshold. Fig. [TJ Pxxi^c) ex- 
hibits four peaks, and so on. By now, such a behav- 
ior (two peaks in Pxxi'-^c) dependence for small enough 
pI^j.) was rep orted in a number of experimental papers 
Refs. la-KLSl. On the other hand, theoretical numeri- 
cal studies^"— aimed at revealing levitation on a mi- 
croscopic level, are less conclusive^^. The only estab- 
lished fact is the tendenci^^— for floating up of delocal- 
ized states upon decreasing Wc in the strong-field domain, 
where Landau levels are still well-defined. This tendency 
is due to disorder-induced mixing of the neighboring well- 
resolved Landau levels. 



B. Network models 

First numerical verification^ of one-parameter scaling 
Eq. ID) was performed within the Anderson model^. 
Physically, this model corresponds to realizations of dis- 
order in which scatterers have random strength, while 
positional disorder is eliminated by placing scatterers on 
the lattice. In other words, the phase acquired by elec- 
tron between two subsequent scattering acts is assumed 
to be the same. Important is that the minimal model^, 
in which disorder is characterized by a single dimension- 
less parameter, spread of the site energies in the units 
of bandwidth, captures all features of a generic random 
potential. 

Another minimal description of disorder is at the 
core of scattering approach to localization introduced 
by B. Shapiro'^-. Microscopic realization of the network 
Ref. t36i requires restricting the electron motion by a con- 
fining potential, as illustrated in Fig. [5] This confining 
potential ensures that there are only four possible out- 
comes of scattering, which takes place at the nodes. In 



contrast to Anderson model, all scatterers at the nodes 
of network are assumed identical, while positional disor- 
der is maximally strong. This is achieved by assuming 
that phases, accumulated between the neighboring nodes, 
are completely random. Within the network model de- 
scription, a physical parameter, (fcj^Z)"^, is emulated by 
1 — T, where T is the transmission of the node. One 
of the apparent successes of the network model descrip- 
tion of disordered systems was the demonstration^ of 
the zero-field localization-delocalization transition in 2D 
system with spin-orbit scattering. Formal origin of this 
transition is the change of sign of P{g) in Eq. ([T|) in the 
presence of spin-orbit scattering^ 

Network-model approach is especially well suited for 
the description of the quantum Hall transition in a strong 
magnetic field. As was pointed out by Chalker and 
Coddington^ (CC), in this case, unlike zero magnetic 
field, the links of the network acquire a natural phys- 
ical meaning, namely, they coincide with equipotential 
lines^"— of the hare smooth (on the scale of magnetic 
length) random potential. This is due to the field-induced 
quenching of kinetic energy of electron, rather than due 
to artificially imposed confining potential. In addition, in 
strong magnetic field, the motion along each link, repre- 
senting the drift of the Larmour circle, is unidirectional. 
Nodes in the of the CC network also have a transparent 
meaning: they represent saddle points of the random po- 
tential, where equipotentials come as close as magnetic 
length. Various aspects of the Quantum Hall transition, 
relevant to experimenti^Ti^, e.g., divergence of the local- 
ization radius (scalin g^^i^^'^^ ) , critical statistics of energy 
levels^, mesoscopic conductance fluctuations^"^, point- 
contact conductance?-, were studied theoretically using 
the CC models. 

Testing the levitation scenario microscopically requires 
to construct a minimal weakly chiral network model, 
which captures the physics encoded in the system Eq. ([5]) , 
namely competition between interference-induced local- 
ization and orbital-induced curving asymmetry in the 
scattering to the "left" and to the "right" . Construction 
of such a network and study of its localization properties 
is the objective of the present paper. On the physical 
grounds, the desired description should contain only two 
parameters, kpl and uJcT. Short communication on the 
results reported below can be found in Ref. [ssl . 
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FIG. 3: (Color online) Construction of weakly chiral network 
model, a: Node scattering matrix is decomposed into four 
"reflectors" and a "junction", b: Final network model upon 
combining two reflectors on a given link into one reflector. 



scattering at the node is described by 4 x 4 unitary scat- 
tering matrix. This matrix can be parameterized by three 
independent numbers, e.g., as follows: 
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where reflection, r, transmission, t, and deflection coeffi- 
cients, di, ^2, are all real and satisfy the flux conservation 
condition 



r"^ + t'^ + dl + dl = 1. 
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Isotropy requires that d ^ di = d2, so that the scatterer 
is characterized by only two independent parameters, say 
t and d. The model belongs to the orthogonal symmetry 
class if the phase, ipi2, accumulated upon propagation 
1 — ?> 2 between the scatterers 1 and 2, Fig. [21 is equal 
to the phase, Lp2i, accumulated upon propagation 2 ^> 1. 
Otherwise, it belongs to the unitary class. 

Full localization of the eigenstates in 2D, predicted by 
the scaling theory^, manifests itself in vanishing trans- 
mission of the network for all sets of t and d. The degree 
of localization is governed by the conductance^! 
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II. REFORMULATION OF NON-CHIRAL 
NETWORK MODEL 

We achieve the goal of constructing a minimal weakly- 
chiral network model in two steps. First we reformulate 
the standard non-chiral network model Fig. [2]by separat- 
ing each node into the regions with backscattering and 
left-right scattering. As a second step we incorporate 
weak chirality in the form of imbalance between scatter- 
ing to the left and scattering to the right. 

In a standard non-chiral network model Fig. [5] the 



lid and t are small, localization is strong (within one pla- 
quette). Conversely, for t close to 1, electron changes the 
direction of propagation after ~ — ^ 1 scattering 
acts, which corresponds to high conductance, kpl ^ 1, 
and exponentially large localization radius. It is impor- 
tant to mention one particular case, namely, t <C 1 and 
d ~ 1/^/2. According to Eq. (fT2|) . this case of almost 
complete deflection, with weak reflection and transmis- 
sion, should correspond to strong localization. This is, 
however, not the case. In fact, the states become pro- 
gressively delocalized as d approaches the value 1/^/2. 
This conclusion can be drawn from Ref. [591, where the 
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corresponding limit of the network model has been stud- 
ied. Since the case d « l/-\/2 will play an important 
role in our construction later on, we discuss a seeming 
contradiction to the scaling in this case in Appendix. 

Note in passing, that the parametrization Eq. (jlOp re- 
flects the Born scattering (for di = ^2), where the prob- 
abilities of deflection, forward, and backward scattering 
are independent of the direction of incidence. This is 
certainly not the general case of potential scattering. In 
fact, one can relax the condition di = d2 for a given node 
and ensure global isotropy by requiring that the scatter- 
ing to the left and to the right are equally probable on 
average (over the nodes). 

Our reformulation of the network model is illustrated 
in Fig. [3l We start with specifying scattering matrices 
at the nodes of Fig. [5] as systems of four "reflectors" 
and one junction, Fig. [3l Each reflector mixes only two 
channels on the corresponding link, so that its scattering 
matrix 



of the scattering matrix with appropriate symmetry. Let 
us denote with 0; the phase accumulated between the 
junction and reflector on the link i, i = 1,..,4, Fig. 
[5^. For potential scattering with time-reversal symmetry 
preserved, the phase accumulated between the reflector 
and the junction is the same 0i, as for propagation be- 
tween the junction and reflector. Then the effective node 
scattering matrix assumes the form 
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is 2 X 2, where po is the power reflection coefficient. The 
junction, Fig. |4l does not transmit or reflect incoming 
waves, but rather scatters them either to the left or to 
the right with equal probability, 1/2. The corresponding 
4x4 scattering matrix has the form 
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It is important to demonstrate that the combination of 
a junction and four reflectors which is characterized by 
a single parameter, p^, corresponds to the effective node 



where the deflection, transmission, and reflection coeffi- 
cients are defined as 



D,^^ = ^^^^ + poe'(2^-^--*^-)) , (16) 



r,,=2z^^°ii^sin(^.-0,)e-^ 
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+2 + 2poe^'^] . (18) 




a b 

FIG. 4: (Color online) a: Microscopic realization of the junc- 
tion matrix Eq. p4|) . The junction is defined by the con- 
finement potential (shaded area); rhomboidal scatterer at the 
center ensures that each incident wave is deflected only to the 
left or to the right, b: Scattering on a link, described by 
matrix Eq. (|2ip . can be modelled with a point contact. 



In Eqs. (fT6)) - (|T8| . ip = J2't'i is the net phase; normaliza- 
tion factor is defined as 



Af = 2 + 2p2e2i'/' + po(e^"^' + e^''l'')ie^''f'^ + e'^"*'*). (19) 



One can check that the matrix Eq. (|15p is unitary. 
For arbitrary phases, (pi, it does not reduce to the form 
Eq. (|10p with equal reflection, transmission, and deflec- 
tion probabilities for all directions of incident channels. 
Global isotropy is restored upon averaging over (pi. For 
example, if we choose a particular set, (pi = (j)2 = 7r/4, 
(pz — <Pi = — ^/4, the matrix Eq. ([15]) assumes the form 
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This choice favors deflection \d2 \ over deflection The 
asymmetry between \d2\ and is compensated by re- 
aUzation in which (pi = 4'2 — —t^/^, 4'3 — 4'i — '^1'^] 
for this reahzation, |c?i| and 1^2! switch places. It is seen 
from Eq. (IT5|) that presence of reflectors eliminates the 
singular character of junction matrix Sq by restoring fl- 
nite forward and backward scattering probabilities. In 
other words, the use of the matrix Eq. instead of 

node matrix Eq. (flUl) would reveal full localization of 
electron states at any po- However, for small poi the 
parametric space of matrices Eq. ([15]) is restricted to 



small transmission and reflection, \t\ 



Po, the 



domain where the mean free path Eq. (1121) is ~ 1. De- 
spite what Eq. ((T2|) predicts, we will get large values of 
localization length, ^{po), in the domain pq <^ 1. 

We complete our reformulation of the non-chiral net- 
work model by observing that the random phases on the 
links between two neighboring reflectors can be incorpo- 
rated into This allows one to combine the two reflec- 
tors on the same link into a single effective 2x2 scatterer 
on this link, as it is shown in Fig. |3}d. The corresponding 
effective scattering matrix, 
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has the same form as Eq. (IT3l) with 
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The resulting network consisting of junctions at nodes 
and effective scatterers on the links is shown in Fig. |3Jd. 
On the microscopic level, the node Sq corresponds to the 
junction with conflnement shown in Fig. |4^, while the 
link matrix V corresponds to point-contact confinement. 
Fig. Hd. 



III. ORBITAL ACTION OF A WEAK 
MAGNETIC FIELD 

As was discussed in the Introduction, weak magnetic 
field in which delocalization transition takes place, is 
already strong enough to drive random phases on the 
links of the network Fig. [3]d into the unitary class. We 



also need to incorporate the Lorentz-force effect of mag- 
netic field. For free electrons, the Lorentz force curves 
their trajectories. In the network Fig. ^jp it affects the 
properties of junctions only, leading to imbalance be- 
tween deflection to the left and deflection to the right. 
Note, that general properties of a four-terminal junction 
in magnetic field were previously studied in Refs. l60l - 
I62I for various forms of confinement potential in relation 
to experiments^^— on the Hall quantization in narrow 
channels. 

In order to incorporate orbital action of magnetic field 
into the general network Fig. [21 one has to place a weakly 
chiral S'-matrix into each node. A possible form of such 
5-matrix is 
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where the complex coefficients have absolute values 
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If magnetic field whirls electrons, say, to the right, one 
has D > d. For example, one can choose the following 
realizations of the above chiral matrix 
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provided that r, d, t, and b± are positive real numbers 
{t > r) with 



t - r 
t + r' 
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We see that the difference between the scattering proba- 
bilities to the right and to the left is \b+\^ — l^-P- This 
difference is non-zero, as follows from first identity in Eq. 
([^51) . Note that for the particular choice Eq. ([25)) . the 
ratio, r/t, "controls" the magnetic field strength. Indeed, 
for r = we have 6+ = 6_. 
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In general, dimensionless Hall resistance of a junction 
is expressed via field- dependent elements of the matrix 
Eq. (Uni) as follows^-^ 



Rh 



2{dl - dl) 



(2<2 



dl 
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Note that for the weakly chiral matrix Eq. (|25l) Rh is 
the same for all directions of incidence. Naturally, the 
degree of bending action of magnetic field is represented 
by the difference, — |<i2p)- The advantage of our 

reformulation of the network model, described in Sec. II, 
is that controlled chirality can be incorporated into the 
network in a natural way upon replacement of the non- 
chiral matrix iSq, Eq. (IT4l) . by 
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The matrix S is parameterized by a single number, q ~ 
dl, which varies between q — Q and g = 1. The Hall 
resistance is expressed via q as follows: 



Rh — 

It is an odd function of 
2 
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i.e., the difference (1/2 — q) can be viewed as a quantita- 
tive measure of the magnetic field strength. Concluding 
this section, the link matrix V, Eq. (|211) . parameterized 
by single parameter, p, and the node matrices 5, Eq. ()28|) . 
parameterized by single parameter, fully define a min- 
imal network model, Fig. [SJj. The advantage of this 
network is that the disorder and the magnetic field can 
be "tuned" independently by changing p and g, respec- 
tively. We will call this model a p-q model. Localization 
properties of this model are studied below. 



IV. P-Q MODEL: LIMIT OF STRONG 
DISORDER 



In order to get a qualitative insight into the phase di- 
agram of the quantum p-q model, we start with artificial 
limit of strong disorder. To define the strong disorder, 
note that in the original p-q model the values of p and 
q are the same for all junctions and point contacts. In 
other words, distribution functions of the parameters p 
and q are 



f{Pi) = 5{p^ - p), f{qj) ^ S{qj - q). 



(31) 



By a strong disorder we mean the following distribution 
of p and g, 

f{Pi)=pSil-p,) + {l-p)Sip,), (32) 
f{qj)=qS{l-q,) + {l-q)S{q,), 
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FIG. 5: (Color online) Left: In two-channel p-q model, elec- 
tron motion is restricted to the spaces between forbidden re- 
gions An,m- The centers of forbidden regions form a square 
lattice. Point contacts on the links describe the backscat- 
tering by disorder; bend-junctions at the nodes describe the 
orbital action of magnetic field. Right: Scattering matrices of 
the junction and of the point contact. 



so that scattering by the point contact and defiection at 
the junction are still p and q on average. However, un- 
like Eq. (PT|) . the point contact reflects fully in p percent 
of cases, and transmits fully in the rest (1 — p) percent 
of cases. Similarly, according to Eq. ([5^. the junction 
deflects only to the right in q^ percent of cases, deflects 
only to the left in (1 — q? percent of cases; in the remain- 
ing 2q{l — q) percent of cases the deflection takes place 
both to the left and to the right depending on incoming 
channel. 

In considering a strong disorder, our motivation stems 
from the original CC model, in which the nodes of the 
network are chiral saddle points. If we introduce disorder 
in the transmission of saddle points similar to Eq. p2p . 
then, at critical energy, 50% of saddle points will fully 
transmit, and 50% will fully reflect. For such a dis- 
order, quantum-mechanical interference becomes irrele- 
vant. However, the state will remain critical®^, separat- 
ing the phases with Uxy differing by 1. In this limit of 
strong disorder, quantum delocalization is replaced by 
classical percolation transition occurring at the same en- 
ergy. Our expectation, which will be later supported by 
numerical simulations, is that similar to CC model, con- 
sidering the limit of strong disorder of the p-q model will 
yield the positions of classical delocalization transitions, 
which coincide with the positions of quantum delocalized 
states in the original p-q model. 

The realization with Pi = \ corresponds to a "closed" 
point contact, which reflects incoming waves from both 
directions. Classically, presence of such a reflecting bar- 
rier can be interpreted as a bond installed between the 
neighboring forbidden regions, An^m and A„+i_,„, of con- 
fining potential. Fig. [S] Below we will refer to this bond 
as a p- bond. Similarly, we introduce q- bonds, installed 
between the forbidden regions, An^m and An±i^m±i in 
Fig. [51 Then, deflection only to the right corresponds 
to two crossed q- bonds, deflection only to the left cor- 
responds to the situation when all four forbidden regions 
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FIG. 6: (Color online) Limit of strong disorder. The centers of 
forbidden regions, A„,m and An,m-i, are connected by the p- 
bond, while the centers of forbidden regions, A„-i^m-i and 
An,m, are connected by a q- bond. Electron delocalization 
reduces to the percolation on the lattice consisting of p- and 
q- bonds. 



■^n,m, defining the junction, are disconnected. One right- 
diagonal q- bond describes the situation when the junc- 
tion deflects the fluxes incident from the left and from 
the right channels to the left, and fluxes incident from 
the up and down to the right. Similarly, one left-diagonal 
q- bond signifies reflection from the left and right chan- 
nels to the right, while the up and down channels are 
deflected by the junction to the left. Thus we arrive at 
the percolation problem of connectivity of the forbidden 
regions An^m via p- and q- bonds, see Fig. [51 Some limits 
of this problem are transparent. For example, for small 
p and g, the A„_,„ regions are mostly disconnected. Also, 
for p > 1/2, global connectivity exists for any q. It is 
intuitively clear that adding small portion of q- bonds 
facilitates connectivity and shifts the position of percola- 
tion transition fromp — 1/2 to lower values of p. Quanti- 
tatively, we will search for the boundaries of percolation 
transition on the p, q plane by employing the real-space 
renormalization group approach to the 2D percolation^. 




a b 

FIG. 7: (Color online) a: Illustration of the real-space 
renormalization-group procedure. The connectivity of a su- 
perbond defined by five p- bonds is enhanced by q- bonds, 
Eq. p9|) . b: Percolation picture at small p. p-bonds "assist" 
the connectivity of q- bonds, Eq. H41[) . 



A. Real-space renormalization-group analysis of 
the percolation problem 

The original approach of Ref. IH applies when g = 0, 
i.e., when only p- bonds connecting the sites of A„,m 
square lattice are present. Within this approach, five 
p- bonds are replaced by one superbond, as shown with 
full lines in Fig. [7^. Probability, p\ that superbond 
is present, is expressed via probability that p- bond is 
present, as 

P' = + Pi{p) + PM + P2{P) + Pl{p) + P0{P), 

= Up). (33) 

where P„ (p) is the partial probability that the superbond 
is present when n original bonds are present. A simple 
counting of variants yields 

Pf,{p)=p\ Pi = bp\l-p), F3 = 8p3(l-p)2, 

P2{p)^2p'{l-pf, Pi=0, Fo-0. (34) 

Remarkable feature of the transformation Eq. is that 
its fixed point, p' — p, coincides with the exact bond per- 
colation threshold, p = 1/2, on the square lattice. The q- 
bonds, shown with dashed lines in Fig. [T] We will incor- 
porate them into renormalization group transformation 
assuming that their role is the enhancement of connec- 
tivity of the superbond. This enhancement occurs differ- 
ently depending on how many original bonds are present. 
For example, if all five or four bonds are present, the con- 
nectivity of superbond is guaranteed even without any q- 
bonds, so that P5 and P4 are not affected by q- bonds. 
When three p- bonds are present, there are two variants 
when superbond does not connect. Then the probability 
that it connects in the presence of q- bonds is given by 

Psip, q) = Psip) + 2pHl - pf [2(7(1 -q) + q^] . (35) 

The second term is the product of probabilities that 
p- connectivity is absent and that q- bonds restore it. 
The factor in the square brackets accounts the fact that 
restoration can happen by installing one q- bond (two 
variants) as well as two q- bonds (one variant). The ex- 
pression for P2{p,q) has a similar structure, 

P2(P, q) = P2{p) +p\l~pf [2 • 4(2g(l -q) + q^)] . (36) 

Overall, there are ten configurations when two p- bonds 
of the superbond are present. Out of these ten, there 
are two variants when the superbond connects; in the re- 
maining eight variants the superbond does not connect. 
The factor in the square brackets in Eq. (1361) describes 
the probability that in these eight variants q- bonds make 
the superbond connect. Note that the q dependence of 
the second term in Eq. ([55]) is the same as in the case 
of P2{p,q). This reflects the fact that in both cases in- 
stalling either one or two q- bonds restore connectivity. 
If there is only one p- bond, it can be either "vertical" 
(one variant) or horizontal (four variants). In the first 
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FIG. 8: (Color online) The line of percolation transitions in 
the {p, g)-plane is plotted from Eq. (|39p for small q (red line) 
and from Eq. (|41|l for small p (green line). Dashed and dotted 
straight lines are the asymptotes, Eq. (|40|l and Eq. (|44|) with 
a = 1, respectively. Crosses show the positions of quantum 
delocalization transition inferred from quantum simulations 
for five values of "energies" , p. 

case, probability of restoring connectivity is a product of 
the probabilities that it is restored both "to the left" and 
"to the right" from the p- bond. If the p- bond is hor- 
izontal, there are three q- bonds that might participate 
in the restoration of the connectivity. This yields 

P,ip,q)=p{l-p)^[{2q{l-q)+q'f 
+4(g(l-g)2 + 3g2(l-g)+q3)]. (37) 

Finally, Po{p,q) is the probability that the superbond 
connects via q- bonds only. All four q- bonds can par- 
ticipate in restoration, see Fig. [7^. In particular, the 
connectivity can go through the upper middle site, lower 
middle site, or both, resulting in 

Po(p, g) - (1 - pY' [2<z'(1 - + 4g3(l -q)+ q^] (38) 

The net probability that the superbond connects is the 
sum of Pr^ip) + Pi{p) and Eqs. dSS]), namely 

f{p,q) = foip) + 2p3(l _p)2[2q- g2] (39) 
+&p^l-pf[2q^q^] 

+p(l - pf [iq + 8g2 _ + q^] + (l - pf [2q^ - q4] _ 

We determine the line of percolation transitions on the 
p,q- plane upon equating f{p,q) to i. Solution of this 
equation is plotted in Fig. [H At small q the p{q) bound- 
ary is linear, 

Pc(g) = ^-^g. (40) 

In the above procedure we assumed that the effect of 
q- bonds is a correction to the percolation over p- bonds. 



Consider now the opposite limit, where q is close to 1/2 
and p is small, so that the percolation is dominated by 
the q- bonds, while the p- bonds constitute a small cor- 
rection. First we note that for p — 0, the regions An^m 
with n + m even and n + m odd are decoupled. Moreover, 
q = 1/2 corresponds to bond percolation threshold in 
both decoupled sublattices. Now adding small portion of 
p- bonds facilitates percolation for q < 1/2. To describe 
this facilitation quantitatively we turn to Fig. [7)3. This 
figure illustrates that, instead of one missing q- bond, a 
pair of one horizontal and one vertical p- bonds can pro- 
vide the connection, and that there are two such variants. 
Resulting shift of the threshold position is determined by 
the condition 

/(p,g)=9' = <Z + 2(l-g^)p2-i, (41) 

where q' is the renormalized probability that q- bond 
connects. The origin of the factor (1 — g^) in Eq. (pij) 
is the following. If the vertical q- bond is missing [with 
probability (1 — q)], the two p- bonds which restore con- 
nectivity can either share a common site (with probabil- 
ity 2p^), or two p- bonds can be connected to each other 
via a horizontal q- bond (with the probability 2p^q). The 
sum of probabilities of these two realizations should be 
multiplied by (1 — q). The resulting from Eq. ()41[) p(q) 
dependence, 

Pc{q)^\l^^, (42) 

is plotted in Fig. [H In fact, as seen from Fig. \8\ the 
asymptotic behaviors found from Eqs. (|39|) and (|4T|) 
match very closely near q — j. 

Eq. (PT|) describes the situation when percolation oc- 
curs in one, e.g., n + m even sublattice, while the n -\- m 
odd sublattice paly an auxiliary role. This picture is 
violated near the degeneracy point p = 0, q — 1/2. 
Upon approaching to this point, both sublattices should 
be treated on the equal footing. Moreover, the position 
of the boundary on the p, q plane is governed not by a 
local arrangement of the bonds but rather by large-scale 
behavior of the clusters in both sublattices. We will dis- 
cuss percolation transitions in this region in a separate 
subsection below, after we establish the relation between 
percolation and electron trajectories in magnetic field. 



B. Implication of the percolation transition for 
transport 

The prime question is: to what extent the connectiv- 
ity of Aij regions studied in the previous subsection gov- 
erns the transport over the regions located between the 
forbidden regions Aij. To emphasize that this question 
is non-trivial, note that, while p- and q- bonds facilitate 
the connectivity, a combination of two p-bonds at a given 
junction and one q- bond through the same junction can 
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FIG. 9: (Color online) a: Structure of electron trajectories 
at p = g = 0. Trajectories either encircle forbidden regions, 
An^rn, (browu liues) , or constitute a chiral edge state due to 
reflection from the boundaries (red line), b: Installing one p- 
(vertical) or q- (diagonal) bond creates a trajectory encircling 
two forbidden regions counter-clockwise, c: Example of a 
finite cluster of p- and q- bonds. Formation of trajectory 
encircling the hull. 



give rise to localized electron trajectory, circling around 
the junction. 

To establish the connection between transport and per- 
colation, we start with the simplest case p — 0, q — 0. It 
is apparent from Fig. [5] (see also Fig. [3^) that in this case 
electron trajectories are closed counter-clockwise loops 




FIG. 10: (Color online) Phase diagram (red line) of the p- 
q model in the regime of strong disorder. Phases (a) and 
(c): electron predominantly encircles forbidden regions An,m, 
moving counter-clockwise (a) or clockwise (c). Coupling of 
neighboring forbidden regions via p- and q- bonds is weak. 
Phases (b) and (d): electron predominantly encircles nodes, 
moving counter-clockwise (d) or clockwise (b). Coupling of 
nodes due to absent p- bonds is weak. Edge trajectory dis- 
appears upon crossing the boundary (red line), (a) — >■ (b) or 
(c) ^ (d). 



motion occurs along the boundaries of forbidden regions 
Aij, at large scales the hulls define the extent of the mo- 
tion. Now we can identify the point of the percolation 



around forbidden regions Aij. In other words, electron threshold at which Aij get connected into infinite cluster 

with the point when the hull trajectories become infinite 
and connect opposite sides of macroscopic sample. 

Essentially, hulls can provide a connectivity through 
an infinite sample exactly at the percolation threshold. 
Below and above this threshold, hulls either do not pro- 
vide a macroscopic connectivity or establish a trajectory 
extending along the macroscopic edges of a sample, de- 
pending on boundary conditions. Being translated into 
the p-q model transport properties, this means that a 
non-zero diagonal conductivity, cFxx 0, is possible only 
at the percolation transition, which separates two insu- 
lating regimes with axx = 0. 

We now choose the boundaries of the macroscopic p- 
q sample passing through the centers of the end regions 



executes counter-clockwise motion along the perimeters 
of Aij. Now let us switch on a single p- bond, say, be- 
tween An^m and An,m-i, Fig. int). It is easy to see that 
installing this bond creates a closed trajectory encircling 
two regions, joined by the bond, in such a way that the 
forbidden joined region remains on the left. Thus, from 
two trajectories along the perimeters of disconnected re- 
gions we get one trajectory along the perimeter of the 
joined region Fig. ^jp. The same happens upon installing 
a single q- bond, as shown in Fig. ^]p. In fact, this evo- 
lution is general: if a portion of p- and q- bonds connect 
several forbidden regions into a cluster, there appears 
a closed trajectory along outer perimeter of the cluster. 
While moving along this perimeter trajectory, the cluster 
remains on the left. For the transport properties of p-q 
model, the outer perimeter, which signifies the most de- 
localized trajectory existing with the given cluster, plays 
a central role. Note that in the percolation theory such 
a perimeter is called a hull. Thus we see that hulls of the 
bond percolation Fig. ^ correspond to the most impor- 
tant electron trajectories of the p-q model, in the sense, 
that, upon approaching the percolation threshold, the 
hulls of big clusters join into even bigger hulls. Extent 
of the region available for electron motion is determined 
by the size of the typical cluster, which is the localiza- 
tion radius of the classical percolation. While locally the 



Aij , as shown in Fig. ^j^, and assume full reflection at the 
boundaries. It can be seen that for p = 0, g = 0, there is 
a macroscopic trajectory spanning near the edges around 
the sample in the clockwise direction. Therefore, in the 
phase a. Fig. 1101 we have nonzero a^y, while cr^x — 
because there is no trajectory through the sample at p — 
q = and in the vicinity. 

Previous consideration pertains to small enough p and 
q. Let us now move along the boundary q = 0. In the 
absence of q- bonds this corresponds to increasing con- 
nectivity, p, in conventional bond percolation problem. 
It is easy to see that the edge state disappears when we 
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FIG. 11: (Color online) Vicinity of the point p — 0, q — 1/2, 
of the phase diagram Fig. 1101 a: Two q-clusters on the same 
sublattice hybridize upon installing a q- bond, b: A q- cluster 
of n + m odd sublattice (upper) and a q- cluster of n + m even 
sublattice (lower) overlap. A joined trajectory (thin blue line) 
is formed upon installing of a single p- bond. The blowup 
illustrates hybridization of trajectories on a microscopic level, 
c: Same as (b), but with two p- bonds connecting q- clusters. 
Blue lines illustrate that the hull trajectory and "internal" 
trajectory are disconnected. 



pass the percolation threshold p — 1/2. Above this point, 
the electron trajectories are strongly localized around the 
junctions (phase b in Fig. I10[) . and describe the clockwise 
motion, unlike counter-clockwise loops in the phase a. 

Now we fixp in the domain (l—p) <^ 1 and move along 
the q- axis. At this point we note that the p-q model 
possesses a duality q ^ {1 — q). One transparent way to 
reveal this duality is to trace the change in strongly local- 
ized trajectories, phase 6, upon changing q by (1— 9). One 
can see that trajectories remain unchanged, while the di- 
rection changes from counterclockwise to clockwise. Such 
a change could be expected on purely physical grounds, 
since transformation q (1 — q) means that the differ- 
ence, (1/2 — q), which represents magnetic field, changes 
sign. Naturally, reversal of magnetic field results in the 
change of direction in which the An^rn regions are circum- 
vented. Therefore there is one-to-one correspondence be- 
tween the strongly localized phase b and phase d. The 
same argument, change of direction of the rotation as a 
result of g — >■ (1 — g) transformation, suggests that the 
phase c is mirror image of the phase a, with the opposite 
sign of (Txy. 



C. Vicinity of the degeneracy point p — 0, q — 1/2 

In the close vicinity, (1/2 — g) ^ 1, p <C 1, of the 
degeneracy point, one has big clusters of q- bonds in 
(n + m) even and (n + m) odd sublattices, which are 
statistically equivalent. Let us start our consideration 
from some q = qo < 1/2 and p ~ 0. Connectivity via 
q- bonds can be enhanced either by adding p- bonds or 
by shifting q closer to 1/2. We argue that, close enough 
to the degeneracy point, both operations are equivalent. 



This is because the number of "even" q- bonds, i.e., the q- 
bonds from the (n -|- m)-even sublattice, and the number 
of "odd" q- bonds, i.e., the q-bonds from the (n-l-m)-odd 
sublattice, involved in a typical cluster, is the same. Once 
this "equal participation" of two sublattices is achieved, 
it is preserved upon increasing both p and q. The value 
of p necessary to achieve this equal-participation regime 
obviously depends on proximity of q to 1/2. A way to 
estimate this necessary p is based on the following rea- 
soning. The spatial separation between neighboring p- 
bonds is ^ 1/VP- If this separation is smaller than the 
typical size of q- cluster, ^ |l/2 — gl"'*^'^, on a given sub- 
lattice, i.e.. 



P> 1 



8/3 



(43) 



then p- bonds connect clusters from different subnet- 
works. The way in which equal-participation regime 
sets in is illustrated in Fig. [11] In Fig. [TTk the con- 
dition Eq. P5|) is not met; clusters grow independently 
in each sublattice as q increases. Figs. [TTb.c illustrate 
that a typical clusters on one sublattices is overlapped 
by some other typical cluster from another sublattice. 
Therefore, few p- bonds per cluster, see Eq. (|43l) . are 
sufficient for formation of a unified cluster. Percola- 
tion boundary, pdq), hes above the boundary of equal- 
participation regime Eq. (|43|) . It is reasonable to assume 
that the percolation threshold corresponds to certain por- 
tion of p- and q- bonds per site. Having in mind that in 
equal-participation regime p- and q- bonds are equivalent 
within a factor, the above assumption leads us to 



Pciq) 



1 

2'^ 



(44) 



with numerical coefficient, a ^ 1. This linear behavior, 
as well as the forms Eqs. (|40| and (|42|) . are in agree- 
ment with the condition Eq. (|331), which in fact should 
hold in the whole domain where we expect delocalization 
transitions with the participation of p- bonds. The above 
reasoning, leading to the linear boundary Eq. (j44p is by 
no means rigorous. We were able to come up with more 
compelling reasoning that for the case when half of q- 
bonds are replaced by 1 — g, the percolation transition 
at small p indeed occurs at p ^ (1/2 — g). Note how- 
ever, that one cannot claim that this auxiliary problem 
provides evidence for linear boundary in our p-q model. 
This is because approaching the point (p, g) = (0,1/2), 
can depend of the direction of the approach, so that the 
approach along a path with " balanced" g and 1 — g bond 
can yield a different boundary. 

Moreover, this linear dependence is a source of an ap- 
parent doubling of the critical exponent when one ap- 
proaches the delocalization point along the horizontal line 
p = const ^ 1. 
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D. Doubling of the critical exponent 

The shape of the percolation transition line Eq. p4)) 
for which the condition Eq. P5)) is met, allows one to 
make quantitative predictions about the behavior of the 
localization radius, ^{p, q). The reasoning goes as follows. 
Eq. ensures that the typical clusters on the two 

sublattices with sizes ^(q) are connected by p- bonds. 
Now, if we keep q unchanged and approach the critical 
line, Pciq), from below, the localization radius grows as 



ap,q) 



aq) 



[pc{q) 



,4/3- 



(45) 



While this relation is asymptotically exact near the delo- 
calization line, we assume that it still holds deeper into 
the region (a). Fig. [101 Then, in order to find an explicit 
dependence ^{p,q), Eq. pi)) can be used in Eq. (HSI) . 
yielding 



ap,q) 



e(<z) 



q)-p 



4/3- 



(46) 



Note now that Eq. (I46p can be interpreted as a doubling 
of the critical exponent. As one starts at some point, 
{po,q), such that pq <^ 1 and {1/2 — q) ~ 1 and move 
towards the transition boundary along the horizontal line 
p = Po, localization radius changes as follows 



C(Po,g) 



1 



,<V2 [l~qf'[l^^^^qf 
For (1/2 — g) ^ po/a Eq. (|T7)) assumes the form 



^{po,q) 



1 



8/3 ' 



(47) 



(48) 



which means that the subcritical behavior of ^{po,q) as 
a function of q is similar to the usual \q — gd""*^^, but 
with the doubled critical exponent. Therefore, it is the 
linear behavior Eq. pi)) which leads to the doubling. On 
the other hand, as we have reasoned in the previous sub- 
section, the linearity of the boundary of the percolation 
transition is a consequence of the "interaction" of two 
subnetworks. 

The remaining question is the behavior of ^(p, q) along 
the line q = 1/2. To address this question, note that 
there are two delocalization transitions which take place 
as q is changed along the line p — pq 1; the first one 
a.i q — qci < 1/2 and the second at <? = qc2 > 1/2. By 
virtue of duality, qd and qc2 are related as qci + qc2 = 1- 
For small p, from Eq. (|44l) we have 



qci[p) = 7; — ' 

2 a 



qc2ip) = l + ^. (49) 
2 a 



In the domain g > 1/2 the ^{p,q) dependence can be 
found using the above arguments 



C{pQ,q) 



1 



q>l/2 



l]4/3 
21 



Po 



,4/3- 



(50) 



RecaU that in Eq. (|47]), the factor [1/2 - g]"^/^ appears 
as a length of the unit plaquette in the percolation over 
p-bonds. Then the second factor, [1/2 — po/a — q]^"^^^, 
can be regarded as a localization radius associated with 
the delocalization transition at (7ci(po)- The general ex- 
pression which respects the duality and contains as limits 
Eqs. (gZl and ^ reads 



ap,q) 



1 



£„„|4/3|1 
a y I I 2 



(51) 



In particular, along the line q — 1/2 this expression pre- 
dicts the following behavior of localization radius 



aP) 



1 



oc 

9=1/2 p 



8/3 ■ 



(52) 



E. Consequences of the doubling 

Below we will demonstrate numerically that Eq. ()51|) 
applies also to the quantum delocalization upon replace- 
ment 4/3 by the quantum critical exponent of localization 
radius. Here we would like to emphasize the similarity 
between the quantum version of Eq. (j51|) and energy de- 
pendence of the localization radius in the case of close, 
e.g., spin-split delocalized states in strong magnetic field. 
Behavior of similar to Eq. (|51l) . was conjectured in 
Ref. [zO and demonstrated numerically in Refs. [71II72I . 
In our case, two close delocalized states correspond to 
opposite directions of magnetic field, and thus have op- 
posite chiralities. On the contrary, in Refs. [70l - l73i chi- 
ralities of both delocalized states, which are close in en- 
ergy, are the same. More microscopic demonstration of 
doubling in the system with two close in energy delocal- 
ized states with the same chirality can be found in Ref. 
Pt^ . The situation considered in Ref. was two layers 
with smooth random potential coupled by tunneling with 
amplitude, to. For a single layer, the transmission of a 
saddle point is given by T{E/r) = [H-exp(i;/r)]-i. p^^. 
two layers, tunneling allows to bypass saddle points. In- 
stead, the role of a saddle point is played by the region 
where equipotentials from different layers come close to 
each other, but do not intersect. It was demonstrated 
in Ref. uA that for E outside the interval (— ioj^o) the 
transmission of such an effective saddle point is given 
by f{E/t) = T(£;Vf2), where f depends on to as 

1/4 



t^'* . Therefor, while in the vicinity of delocalized states 
E = ±to the critical behavior of localization length is 
^ oc |£' ± to\~'^, outside the interval {—to, to) we have 
^ oc {E'^)~'^ , which mimics the doubling of the critical 
exponent. 



F. Physical interpretation of the phase diagram 
Fig. [H 

In physical terms, the p-q boundary in Fig. [TU] relates 
the backscattering probability, p, and magnetic filed. 
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FIG. 12: Illustration of the transfer matrix, T, Eq. (|53[l . of 
the p-q model. Two q- bonds correspond to the first matrix 
in the product Eq. (|53|) : two p- bonds correspond to the third 
block-diagonal matrix in the product. 

(1/2 — g), at which delocalization transition takes place. 
The parameter p can be viewed as a measure of disor- 
der and also as a measure of electron energy^ Ep] ob- 
viously, p decreases monotonously with increasing Ep. 
Then, the domain q < 1/2 of the phase boundary in 
Fig. [TU] translates into the low-field region, lUcT < 1, 
of the dependence Eo{ujc) in Fig. [TJ Correspondingly, 
the domain q > 1/2 maps onto Eq{uJc) dependence with 
reversed sign of the magnetic field. More detailed cor- 
respondence between Figs. [T] and [TU] can be established 
upon identification of (r^y values for different phases. As 
we demonstrated above using percolation language, there 
is no edge state in the central phase, which includes re- 
gions b and d in Fig. 1101 Since this phase represents 
the region Ep < Eo{ujc) in Fig. [U this phase should 
be identified with the Anderson insulator. We have also 
demonstrated, see Fig. [3^, that there is one chiral edge 
state in the phase a of the phase diagram. Thus, in the 



region Ep > Eq{wc) in Fig. [T]we have axy = 1, so this 
phase is a quantum Hall insulator. In the region p > 1/2 
there are no delocalized states in Fig. [TOl This region 
translates into the strongly localized regime, Ep < 1/r, 
i.e., below the minimum of the _Bo(wc) curve. 

Linearity of the Pc{q) boundary can be rewritten in 
terms of observables. Upon identifying p with l/{kpl) 
and (1/2 — q) with uJcT, the position of the boundary can 
be presented as {kFl){oJcT) — const. This quantifies the 
levitation rate, uJcT oc l/i^kpl)- Curiously, it coincides 
with the prediction of scaling theory, Eq. ([5]). 

For conclusive confirmation of the levitation scenario it 
should be demonstrated that the p-q phase boundary re- 
tains its shape in the presence of quantum interference. 
In fact, the quantum and "percolation" boundaries al- 
most coincide. The evidence for that will be presented in 
the next Section. 

V. NUMERICAL RESULTS FOR QUANTUM 
DELOCALIZATION 

A. Transfer matrix 

Scattering matrices on the links and at the nodes are 
given by Eqs. (PT|) and (|^ . respectively. With regard to 
numerical simulations, the p-q model is quite similar to 
the models with mixing of two co-propagating channels 
on the links, studied in Refs. 28,29,71J2. 

The transfer matrix, T, at each node of the network is 
a 4 X 4 matrix, which transforms four amplitudes on the 
left into four amplitudes on the right. We incorporate 
the link p-matrices and junction q-matrices into T in the 
way illustrated in Fig. [T^l It follows from Fig. [T^ that 
T matrix can be parameterized as 
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(53) 
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The second and the fourth matrices in the product de- 
scribe random phases acquired between the point contact 
and the junction. 

Once the transfer matrix is specified, one can find the 
transfer matrix of a slice with AI nodes in transverse 
direction. As usual, to study the critical properties, nu- 
merical simulations are performed on a system with fixed 
width AI — 2*^, then repeated for 4,5, 6, 7 and sometimes 
even k = 8. By multiplying transfer matrices for a stripe 
with N slices and diagonalizing the resulting total trans- 
fer matrix, it is possible to extract the smallest positive 
Lyapunov exponent Xm/2 (the eigenvalues of the trans- 
fer matrix are exp(AiA^)). The localization length, ^m, 



is proportional to 1/\m/2- The state is identified as a 
critical one when the renormalized localization length, 
£,m/M, becomes independent of the system width, AI. 
We consider this as a criterion for a transition point, 
that determines the value q = qc, for a fixed value of p. 
Practically, q = qc emerges as a maximum in (,m/M vs. 
q for large system widths, when finite-size errors become 
small. 

Close to g = qc, the ratios £,m/AI should satisfy a one- 
parameter scaling 
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which is commonly used to infer the locahzation 
length, ^. 

We first check that our numerical data support 
Eq. ([Ml) : to do so we fit all the data points onto one 
curve according to Eq. ((54)) . The fit is carried out with 
the help of a special optimization program. This op- 
timization program runs different critical g-values and 
critical exponents v and chooses the optimal sets, qc and 
i^, which provides the best agreement with Eq. (|54|) . 

We now briefly describe the optimization procedure. 
The routine determines least-squares polynomial approx- 
imation by minimizing the sum of squares of the devia- 
tions of the data points from the corresponding values 
of a polynomial. The argument of the function fitted by 
the Chebyshev polynomials is Af jq — gd"- To choose the 
optimal pair, q^. and we run this routine for the wide 
range of values q^ and v. For most of the points on the 
critical line the values of (jc, obtained by two methods: 
(?) searching for q at which £,m/M is constant and (m) 
using optimization procedure, agree with each other. 

Note however, that in two limiting cases, p — > and 
p — 1/2, where the data strongly fluctuate, apparent 
discrepancies arise. For very small values of p the off- 
diagonal terms in the transfer matrix are close to (they 
are ^ y/p), leading to strong fluctuations in numerical 
results. Physically, enhancement of fluctuations near p = 

is a result of proximity to two critical points, qc and 

1 — qc, see Fig. [101 where the doubling of critical exponent 
takes place. On the other hand, when p is close to 1/2, we 
have qc close to 0. Strong fluctuations in the data in this 
case is a consequence of small denominators, Xj^fq^ in 
the transfer matrix Eq. (|53l) . Altogether, both methods 
yield close values of <?c(p)- 

The last remark on simulation procedure is on the 
boundary conditions in transverse direction. In CC 
model, the periodic boundary conditions in transverse 
direction are insured upon imposing requirement on the 
structure of transfer matrix of even slices only. The same 
is true for the p-q model. Unlike the CC model, where 
the reflection and transmission at the nodes alternate be- 
tween subsequent slices, the elementary transfer matrices 
in the p-q model, Eq. (|53p. are the same for even and 
odd slices. This equivalence is a result of symmetry of 
a single node with respect to 90° rotations in the p-q 
model. 



B. Zero magnetic field 

Zero magnetic field corresponds to the line g = 1/2 on 
the p, q plane. Above we identified this line with the ver- 
tical energy axis in Fig. [TJ Since the key ingredient of the 
levitation scenario is that all the states on this axis are 
localized, we start with studying localization properties 
along the line g = 1/2. 

To obtain a scahng plot, Ca//M = f{i{p)/M), 
we analyzed 40 data points: four M-values, 
M = 16, 32, 64, 128, and ten p-values, p = 
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FIG. 13; (Color online) A fit of S^m /M data points for the p-q 
model a± q = 1/2 to a one-parameter scaling form: In {(,m/M) 
vs. ln(5/M). 



0.05, 0.1, 0.15, 0.2, 0.25, 0.35, 0.45, 0.53, 0.6, 0.8. By 
plotting ln{^M /M) vs. — InM and shifting points for 
different p to fall on the same line, we get the result 
shown in Fig. [T31 It is seen that the quality of scaling 
is high. The scaling function, /, found from the data 
in Fig. [T31 is shown in Fig. [TH The dependence, 
^(p) inferred in this way is plotted in Fig. [TS]with the 
black line. Our result confirms the expectation that 
the localization length in zero magnetic field increases 
rapidly even in log-scale as p goes to zero. Scaling theory 
predicts the dependence, ln^„ ^ (kpl)^, Eq. ([3]). Since 
we have earlier identified kpl with 1/p, we expect the 
dependence, In^(p) ^ 1/p^. The black curve in Fig. [15] 
falls off with p slower, and can be well approximated with 
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FIG. 14: (Color online) One parameter scaling function, 
^m/M = f{^/M), is plotted from Fig. [T31 
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FIG. 16; (Color online) A fit of £,m/M data points for the p-q 
model at g = 1/2 with time reversal symmetry on the links 
to a one-parameter scaling form £,m/M = f{^/M). 



FIG. 15: (Color online) Localization radius of the p-q model 
in a zero magnetic field (g = 1/2) is plotted vs. inverse "en- 
ergy", (p) for the cases: no spread in scattering strengths on 
the links, p = (p), and no TRS on the links (black), p = (p), 
with TRS on the links (red). Blue and green curves show In^ 
with and without TRS, respectively, plotted for the strong 
disorder in scattering strengths on the links: p randomly as- 
sumes the values 0.01 and 0.99 while the inverse "energy" is 
the average value, (p). Grey curve [ln(^/5o) ~ — f Inp] is the 
result of the percolation treatment Eq. (|52p . The constant 
is chosen to match the data at (p) — 0.8. 



In^(p) = 0.99 — 4 Inp. One possibility to account for this 
discrepancy is that asymptotic behavior is achieved 
at p smaller that 0.05- the minimal p we studied. Note 
however that the total range of change of ^(p) in the 
domain we studied is huge: ^(0.05)/^(0.8) ~ 5.5 x 10^. 

Another prediction of the scaling theory is that the 
presence of time reversal symmetry (TRS) the growth of 
localization radius with kpl is much slower, In^o kpl, 
Eq. ([2]). In the above simulations we assumed that 
magnetic field is zero in the "orbital" sense {q — 1/2), 
but phases on the links were " unitary" . This is because 
the orthogonal-unitary crossover takes place at exponen- 
tially small {1/2 — q). However, in order to relate closer 
our calculation to the scaling theory, we ran simulations 
with TRS on the links restored. This amounts to set- 
ting ip2 = -"^1, ifii = -m, 4>2 = -01, and 04 = -03 
in Eq. (|53| . The scaling function obtained with TRS is 
shown in Fig. 1161 and the corresponding ^(p) is plotted 
in Fig. [13] with the red curve. As could be expected, in 
the strongly localized domain, p > 0.6, there is no dif- 
ference between the unitary and orthogonal cases. Upon 
decreasing p, orthogonal ^(p) indeed grows much slower 
than unitary ^(p), so that ^(0.05)/C(0.8) w 48.5. The 
fact that orthogonal In^(p) extrapolates at p — > to a 



finite value also suggests that diverging behavior sets in 
at p smaller than 0.05. 

In the above analytical treatment of the p-q model we 
considered the limit of strong disorder, which is a strong 
spread in the local values of p with average p — (p) 
fixed. Consideration was based on the understanding 
that strong spread in p eliminates completely the inter- 
ference effects, and thus reduces the analysis of the p-q 
model to the percolation problem, which predicts much 
smaller localization radius, £,{p) ^ 1/p^^^, Eq. ([5^ . In 
order to test this expectation, we incorporated a strong 
spread in p into transfer-matrix calculation. Namely, 
for a fixed (p) we randomly set the values p — 0.01 or 
p = 0.99 on each link. Our results on ^(p) with and 
without TRS are shown in Fig. [12] with the blue and 
green lines, respectively. We see that our expectation is 
confirmed within the region p > 0.5. In this region, the 
two curves with disorder are not sensitive to universality 
class and are well below the curves without disorder in 
p. Moreover, their behavior in this region is in accord 
with prediction of percolation theory. Indeed, percola- 
tion theory predicts C(0.53)/C(0.8) = (0.8/0.53)^/3 ^ 3_ 
From the blue and green curves we get the close values, 
3 and 3.67, respectively. Fig. [TSl also illustrates that for 
p < 0.5, quantum mechanics "wins" over disorder: the 
curves with random p merge with curves without disorder 
in p corresponding to their respective symmetry classes. 

The above results pose an acute question: whether 
the boundary, pdq), of delocalization transitions, which 
was established within the percolation treatment, and 
extends in the region p < 0.5, is preserved in fully quan- 
tum limit. A related question is whether this boundary 
is sensitive to the universality class. We address these 
questions below. 
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FIG. 17: (Color online) Delocalization points for three differ- 
ent sets of p: non-random p, random p with given (p), and 
non-random p with the TRS. 
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FIG. 18: (Color online) Dependence of the critical exponent v 
onp as inferred from optimization procedure, described in the 
text. As p decreases from 1/2 to 0, u{p) grows approximately 
twice. 



C. The line of delocalization transitions and 
critical exponent 

Our main result, Pc{q) boundary without TRS, is 
shown in Fig. [8] with crosses. Five points obtained for 
p ^ 0.1, 0.15, 0.25, 0.35, and 0.45, essentially fall onto 
a straight line, p = 1/2 — q. Deviation from this line 
takes place near p — 1/2, where the data-points match 
well the results of percolation treatment (red line). This 
confirms our expectation that percolation treatment of 
the p-q model correctly predicts position of the quantum 
delocalization transition. Certainly, the critical exponent 
of delocalization transition is different: v — 1 jZ instead 
of = 4/3 for percolation. For smaller p, deviation from 
the percolation treatment (green line) is notable. On the 
other hand, as we explained above, our percolation esti- 
mate is rather rude at small p. We also argued that the 
true small-p percolation boundary should be linear. This 
linearity would mean that the constant, a, in Eq. (j44p is 
a = 1. ' 

More numerical results for p-q boundary are presented 
in Fig. 1171 For comparison, we reproduced the blue 
crosses from Fig. [5] in Fig. [TT] Red crosses demonstrate 
that the boundary is robust against strong disorder in 
backscattering strength, p. In the same way as for zero 
magnetic field, the disorder in p was incorporated by 
randomly choosing local values of p to be 0.99 or 0.01 
while keeping (p) fixed. Coincidence of the two bound- 
aries indicates that delocalization transition is governed 
exclusively by average p, but does not depend on local 
disorder. 

The fact that in the CC model the position of the de- 
localized state does not depend on the spread in local 
transmission coefficients of the nodes (i.e., the spread in 
saddle point heights), is obvious consequence of duality. 
On the other hand, the p-q model does not possess du- 
ality, so that coincidence of the two boundaries is by no 



means obvious. It can be interpreted as an evidence that, 
even in magnetic field, localization properties of the sys- 
tem are governed by zero-field conductance, ~ 1/ (p) ■ 

Change of the universality class amounts to replace- 
ment 1/axx (unitary) by a constant (orthogonal) in the 
first term of the scaling equation Eq. ([5]). If the scal- 
ing theory applies, this replacement should not affect the 
fixed point, axy = n. In other words, only orbital action 
of magnetic field is sufficient to drive the system into 
quantum Hall insulator state^^. We were able to check 
this prediction within the p-q model. Green crosses in 
Fig. [T7] show the position of p-q boundary with TRS. 
Overall, the boundaries with and without TRS coincide. 
Discrepancy at p = 0.1 is likely due to strong fluctuations 
of the data at small p. 

Together with p-q boundary, optimization procedure 
produced the value of critical exponent, v. The depen- 
dence vIp) is shown in Fig. 1181 Horizontal lines in Fig. 
1181 are drawn to illustrate that simulations indicate ap- 
parent doubling of the critical exponent, discussed for 
classical percolation in the previous section. As follows 
from this discussion, the doubling, revealed by numer- 
ics, is a consequence of two close delocalized states. If 
scaling analysis could be carried our in the immediate 
vicinity of a given delocalized state, it would recover the 
conventional value, v 7/3. Indeed, for p = 0.45, where 
two delocalized states are far apart, optimization yields 
1^(0.45) = 2.3, while for p — 0.1, where the fluctuations 
are strong, we get 1^(0.1) = 4.9, which is even bigger 
than 2 • 7/3. What is remarkable about Fig. [TH]is that 
the apparent growth of p upon decreasing p starts quite 
early, e.g., for p = 0.25 we get ly — 2.9. To make sure 
that optimization does not distort the raw data, we have 
checked the scaling manually and reproduced the largest 
and smallest values of v. 
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FIG. 19: (Color online) Triangular version of the "square" p-q 
model, Fig. [5] Point contacts on the links stand for the same 
2x2 scattering matrix as in Fig. [S] while the bend-junctions 
are described by a 6 x 6 matrix Eq. (|64p . The centers of the 
forbidden regions (shaded triangles) constitute a hexagonal 
lattice. 



VI. TRIANGULAR P-Q MODEL 



1 2 3 4 2N-2 2N-1 2N 
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FIG. 20; (Color online) A slice of triangular network of width 
Af = 8 is shown. Three amplitudes on the links to left of 
green dot are related to three amplitudes to the right of green 
dot via matrix X. The upper and the lower boundaries of the 
slice are connected by dashed-dotted lines manifesting that 
the amplitudes on these boundaries are the same by virtue of 
cyclic boundary conditions. Upon passing the white stripe, 
the vector of the amplitudes {Zi \ is multiplied by the matrix 
Li. Upon passing the blue stripe the vector {Zi} is multiplied 
by L2. 



In the previous consideration, electron motion was re- 
stricted to the channels between forbidden regions, An^m, 
with centers residing on a square lattice. This consider- 
ation led us to the phase diagram Fig. [SI containing the 
line of delocalization transition in the p, q- plane. In the 
present section we will demonstrate that the same shape 
of the transition line emerges when the centers of the for- 
bidden regions constitute a hexagonal lattice, as shown 
in Fig. [T21 This figure illustrates that, similarly to the 
square p-q model, backscattering takes place on the links. 
Fig. [Tni also illustrates that for hexagonal arrangement 
of forbidden regions, a junction corresponds to the point 
where six such regions come close. 

By contrast to the square p-q model, a wave incident 
on a junction can be scattered not into two but rather 
into three directions. Namely, it can proceed forward, or 
get deflected by the angles ±7r/3, Fig. [121 Recall that 
for the square p-q model the 4x4 scattering matrix of 
a junction Eq. (|28p had a simple form, namely, a di- 
rect product of two 2x2 matrices. Correspondingly, the 
6x6 junction scattering matrix in Fig. [TO] is a direct 
sum of two 3x3 matrices [see Eq. (|64| below]. This 
is a consequence of the fact that the two channels on a 
given link are not mixed by the junction. As a result, 
similarly to the square p-q model, upon switching off the 
backscattering on the links, the network breaks into two 
decoupled fully chiral networks. While for the square p-q 
model each chiral network was of CC type, here each chi- 
ral network represents a "triangular" model introduced 
in Ref. This chiral model, illustrated in Fig. [20l 

is called triangular because the nodes are arranged on a 
triangular lattice. For convenience, we briefly review the 
triangular network model Ref. [t^ in the next subsection. 



A. Fully chiral triangular model 

Microscopic picture underlying CC and triangular chi- 
ral models is illustrated in Fig. [5T] Fully chiral motion 
represents drift of the Larmour circle along equipotential 
lines of smooth random potential. In CC model equipo- 
tentials meet pairwise at the saddle points. As the en- 
ergy, e, passes through zero, the geometry of equipoten- 
tials evolves from reflection to transmission, as shown in 
Fig. [2Tk . The scattering matrix describing this evolution 
near e = has the foim^ 



Sccie) 



1 _ 

^/2 " 

1 
V2 



1 
V2 



V2 



(55) 



At £ = the power reflection and transmission coeffi- 
cients are both equal to 1/2. 

In triangular model it is assumed that smooth poten- 
tial has 120° rotational symmetry. As a result, equipo- 
tentials meet in the groups of three. As e is swept through 
zero, they evolve from reflection (shaded regions discon- 
nected) to transmission (shaded regions fully connected) , 
see Fig. [2Tb . Near e = 0, probabilities of scattering to 
the left, to the right and forward, are all finite. Obviously, 
at e = probabilities of the left- and right- scattering are 
equal to each other. In Ref. 1?^ it was demonstrated that 
these probabilities are equal to 4/9. Correspondingly, 
the probability of the forward scattering is 1/9. At small 
but finite e the form of scattering matrix is dictated by 
e — > — £ duality and flux conservation. Up to terms it 
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is given by 



/ 1(1 



1(1 -£) 



• (56) 



(l-£) |(l + e)/ 



Recall that in the limit of strong disorder, when the black 
regions in Fig. [2Tk are either connected or fully dis- 
connected, the CC model reduces to the bond percola- 
tion problem on a square lattice. Correspondingly, the 
strong-disorder limit of the chiral triangular model is the 
site percolation on a triangular lattice^^. In the strong 
disorder limit, e = corresponds to the equal portion 
of present and absent sites. In accord to well-known re- 
sult, this site percolation problem possesses a property 
of self-duality^^. 

In Ref. 1^ the critical exponent of the triangular model 
Fig. 1201 was inferred from the real space renormalization 
group analysis. The result, « 2.3 2.76, agrees with 
simulations of the CC model reported in the literature. 
Here we present the result of numerical simulations of 
the triangular model. Simulations use the slice transfer 
matrix 



n 



(57) 



Operators Li and L2 act in "white" and "blue" stripes 
in Fig. [20I respectively. The operator Li performs 
transformation of the vector of amplitudes, {Z^'-"-'}, into 
while L2 performs transformation of {Z^"^} 
into see Fig. [201 For a particular stripe width. 
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FIG. 21: (Color online) Evolution of the equipotential lines: 
(a) in the CC model, (b) in the chiral triangular model. 



M = 8, the matrix forms of Li and L2 are the following 
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(58) 
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(59) 



with dots standing for zeroes. In the above relations, 
{xij} form a 3 X 3 matrix, X, which is the node transfer 
matrix corresponding to the scattering matrix Eq. (|56p . 
and up to terms has the matrix form 



/2(l + £) -2(1 + e) 1 \ 

2{l + e) -3 2(1 -e) 
V 1 -2(1 -e) 2(l-e)/ 



(60) 



Specific form of L2 accounts for the cyclic boundary con- 
ditions in the vertical direction. Matrices 
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FIG. 22; (Color online) Numerical results for fully chiral tri- 
angular network model. A fit of data to a one-parameter 
scaling form, ^m/M = /(eAf^''") yields the critical exponent, 
v ^ 2.36. 
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FIG. 23: (Color online) Numerical results for the line of de- 
localization transitions in triangular p-q model. The point 
q = Q corresponds to p = pc = 0.6527- threshold of bond 
percolation on hexagonal lattice. 



and 



,(") 



diag{e 



}■ 



(62) 



account for the random phases on the hnks. 

Similarly to CC model, the critical properties of the 
triangular model were inferred from the scaling analysis 
of the Lyapunov exponents. The scaling plot is shown in 
Fig. [22l A high quality scaling was achieved for v = 2.36, 
which is in a good agreement with simulations of the CC 
model. 



a function of energy, e. By contrast, in the p-q model, 
describing low magnetic fields, the scattering matrix at 
the node depends on "magnetic field", 1/2 — q, while the 
energy dependence enters via the backscattering proba- 
bility, p. Still, the structure of the scattering matrices 
of the CC model and the p-q model are the same. It is 
also the case for the triangular p-q model, for which we 
choose the following form of the q- dependent scattering 
matrix 



/ l-q q{q-l) 



S 



g + 1 



q 1-9 9(5'- 1) 

\q{q~l) q l^q J 



(63) 



This matrix is unitary for all < g < 1- It is critical 
aX q — 1/2. Indeed, as follows from Eq. the ratio 

of probabilities of scattering to the left and to the right 
is q/{l — q), so that at q = 1/2 these probabilities are 
equal. Their values are 4/9 in agreement with Eq. (j56p . 
With this parametrization, the 6x6 junction scattering 
matrix in Fig. \W\ acquires the form 



/ l-q q{q-l) q \ 

l-q g 

g l-q q{q-i) 

q l-q q{q-l) 

g(g-l) g l-q 

\q(q-l) q l-q J 



(64) 



B. Numerical results for the triangular p-q model 

In the CC model, which is based on the picture of 
equipotentials, the scattering matrix at a node, Sccj is 



Since backscattering in triangular p-q model. Fig. 1201 
takes place on the links, it is described by the same 2x2 
matrix Eq. pip as in the square p-q model. Thus, by 
analogy to Eq. ()53p . the T- matrix of the triangular p-q 
model has the form 
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(65) 



where the diagonal matrices 

Di =diag{e*^\---,e*'^''} 



(66) 



and 

D2 =diag{e"^S---,e*'^«} (67) 
account for the random phases on the links. The matrix. 
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{xij} = X, has the explicit q- dependence 



9 



X(g) = 



1 \ 



1 q-q+1 1 
9 9(9-1) g-1 



(68) 



and {xij} = X are given by the relation 



(69) 



In the square p-q model, two prominent points on the 
p,q plane were (^,0) and (0, |). They corresponded to 
the bond percolation over p- and q- bonds, respectively. 
Similarly, in the triangular p-q model, the point (0, ^) 
is distinguished. At this point, due to the absence of p- 
bonds, two decoupled q- subnetworks undergo the site 
percolation. In triangular p-q model, the counterpart of 
the p- bond percolation at (5, 0) of the square p-q model 
is a p- bond percolation on a hexagonal lattice. This is 
because the centers of forbidden regions in Fig. [TOlconsti- 
tute a hexagonal lattice. Thus, the second distinguished 
point for triangular p-q model should be {pc, 0), where 



Pc 



1 



2 sin ( — 
V18 



0.6527 



(70) 



is the threshold of the bond percolation on the hon- 
eycomb lattice^. Our numerical simulations based on 
the transfer matrix Eq. (j65p confirm this expectation. 
The points of the delocalization transition, shown in Fig. 
[23l follow the line which smoothly connects the point 
(0.6527,0) and (0, ^). The general shape of the line is 
quite similar to the above results. Fig. El for the square 
p-q model. In particular, in the domain of vanishing 
magnetic fields, (1/2 — q) ^ 1, the transition boundary 
is again linear, in agreement with prediction Eq. ^ of 
the scaling theory. 

In conclusion of this section we would like to make the 
following remark. Simplification of the node structure 
in the square p-q model was achieved by choosing the 
matrix Eq. (I28|) which captures orbital action of mag- 
netic field, but does not allow forward- and backward 
scattering, which, instead, takes place on the links. Sim- 
ilarly, in the triangular p-q model, the junction matrix 
Eq. ([M)) restricts the scattering options for incident elec- 
trons. Namely, out of six possibilities, the electron can 
be scattered only into three channels. Again, similarly 
to the square p-q model, it can access the three other 
channels upon backscattering on the links. 



VII. STRONGLY LOCALIZED REGION: 
IMPLICATIONS FOR INELASTIC TRANSPORT 

The presence of a minimum in the line, En{uJc), of 
the delocalization transitions, Fig. [1] manifests itself in 
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FIG. 24: (Color online) Levitation plot for n = 0; same as 
lower curve in Fig. [T] For electron densities corresponding to 
Ef slightly below the minimum, inelastic magnetoresistance 
has a deep minimum near magnetic field, Uc = i^c- 



specific behavior of low-temperature magnetoresistance, 
Pxx{'-^c)- For high enough electron densities, pxx{'-Oc) is 
temperature independent at two (low and high) distinct 
values of magnetic field. This T- independence is a sig- 
nature of the delocalization transition. Different exper- 
imental groups have studied details of the behavior of 
Pxx{i^c,T) near the high-field transition22r^j low-field 
transition22ii^, or both transitions^"—. For Fermi level 
position below the minimum in Fig. [1] the system re- 
mains localized upon increasing ujc, see Fig. [24l How- 
ever, even in this localized regime, proximity of ujc to the 
position of minimum, w^, (Fig. I24p should manifest itself 
as a precursor of delocalization in the inelastic transport. 
Such a precursor of delocalization has actually been ob- 
served in the early papers^i2£ in the form of a minimum 
in Pxx{'-^c), in the regime of the variable-range hopping. 
This minimum reflects the increase of localization radius, 
^(wc), near uJc = w^- The lower is the temperature, the 
deeper is minimum, as follows from the Mott's law 



In [pxxii^c)] oc 



1 



1/3- 



(71) 



Theoretical studies of variable-range- hopping magnetore- 
sistance pertained to deeply localized regime, where in- 
terference effects in a single hopping act constituted 
a small but singular in lUc correction to the tunneling 
probabilit y^^ 1^^ . The correction is small because the tun- 
neling probability is greatly reduced if electron under- 
barrier trajectory deviates from the straight line. Ac- 
cording to Refs. [91I, m, interference responsible for mag- 
netoresistance occurs between different virtual forward- 
tunneling trajectories. The other early theory Ref. (9^ 
of negative hopping magnetoresistance was based on the 
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FIG. 25: (Color online) Localization radius in the strongly 
localized region p > 1/2 is plotted vs. magnetic field. The 
point q = Q represents the field, ujc = lo'c in Fig. [24l Red and 
black curves for p = 0.53, and blue and green curves for p = 
0.6 corresponding to TRS and no TRS, respectively, merge 
upon increasing magnetic field. Light blue curve corresponds 
to p = 0.8. There is no difference between TRS and no TRS 
for this deeply localized energy. 



following reasoning. Weak magnetic field, by changing 
the universality class and thus inducing delocalization of 
states vi^ith high energies, > l/r, causes some growth of 
^ for the states in the deep tail. Note that both theories, 
as well as later theory Ref. 94, were based on the phase 
rather than orbital action of magnetic field. 

By now there is no theory describing the behavior of 
localization radius, ^{Ep,uic), in the domain of interme- 
diate fields Wc ^ ~ l/r and electron densities ( Fermi 
energies, Ep) at which electron states are localized but 
not strongly, so that EpT ^ kpl ^ 1. On the other 
hand, the p-q model offers a unique possibility to study 
localization length in this domain. To this end, we stud- 
ied numerically ^{p, q) in the region p > 1/2 on the p, q 
plane. The dependence ^(p, q) can be translated into the 
S^{Ep,uJc) dependence. We find that the phase action of 
magnetic field does lead to a certain increase in ^ at low 
fields. However, further increase of Wc, when the orbital 
action sets in, causes a much stronger delocalization ef- 
fect. 

Our numerical results for ^(p, q) are presented in 
Fig. [551 for the values p = 0.53, 0.6, and 0.8. In or- 
der to illuminate the role of the mechanism Refs. |9l|, [92l . 
simulations were performed with and without TRS on 



the links. The values of £,{p,q) were inferred from the 
scaling analysis of /M, similarly to the case of zero 
magnetic field. For p ~ 0.8, the difference between TRS 
and no-TRS is negligibly small. Still the interplay of or- 
bital effect and interference leads to enhancement of f 
from weak, q = 1/2, to strong, q = 0, magnetic fields, 
by a factor of 2. From Eq. ([7T|) it follows that the cor- 
responding drop-off of log-resistance is 2^/^, and \n{p) is 
approximately linear in (1/2 — q). This allows a com- 
parison with the experimental data of Ref. 16|, where gi- 
ant negative hopping magnetoresistance of a degenerate 
2D electron gas was reported. For experimental value 
of ^ = 250 A the low-field log-resistance at temperature, 
T = 0.3K, was 8.7. The net drop-off of log-resistance 
observed was 2.3. This corresponds to the increase of ^ 
by a factor of 1.5. However, at p = 0.8 our numerics sug- 
gests a linear change of ^ with q, while experimentally it 
changes slower. The possible origin of this discrepancy 
lies in the fact that (1/2 — q) is proportional to ujc only 
at low fields. In fact, at higher fields, (1/2 — q) changes 
slower than uJq- 

As p decreases towards 1/2, both phase and orbital 
mechanisms of magnetoresistance become stronger. As 
follows from Fig. the (zero-field) ^ increases from 
TRS to no TRS by a factor 1.47 and 1.82 for p = 0.6 
and p = 0.53, respectively. This increase is a quantita- 
tive measure of the phase mechanism. It should result 
in drop-off of the log-resistance by the factors 1.3 and 
1.5, which corresponds to the resistance drop by several 
times. The orbital mechanism causes a significantly big- 
ger reduction of log-resistance. The p = 0.6 and p — 0.53 
curves in Fig. [^indicate the decrease in ^ between q = 
and q = 0.5 by factors 11 and 60, respectively. Both ef- 
fects are way too strong in comparison with all the ex- 
perimental data on hopping transport reported in the 
literature. This is because p = 0.6 and p — 0.53 cor- 
respond to very large values of localization radius, e.g., 

— 3000A in experimental conditions of Ref. Temper- 
atures required to observe hopping transport with such 
large ^ are unreasonably low. 



VIII. CONCLUSION 

It is expected on general grounds that, as disordered 
system undergoes a quantum Hall transition, its behavior 
is universal, i.e., in the vicinity of transition this behavior 
does not depend on the type of disorder. For quantum 
Hall transition in a strong magnetic field numerics pro- 
vides a compelling evidence for such a universality. For 
example, critical exponent derived from the CC model, 
based on the picture of a smooth disorder, and from the 
simulations^^ for the point-like scatterers are the same. 

Scaling scenario of Ref. HI suggests that regardless of 
the type of the disorder, a single characteristics of the 
system, the Drude conductance a, determines the size of 
the sample, or depending on presence or absence 
of TRS, at which the system becomes an insulator. By 
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virtue of scaling scenario, electron moving on the links 
of the network and scattered at the nodes gets localized, 
due to interference effects, in the same way as realistic 
electron moving on a plane and scattered by random im- 
purities. Correspondence between two systems is estab- 
lished by Eq. (fl^ . 

Pruisken's theory suggests that, in magnetic field, only 
two characteristics, the components axx and a^y of the 
Drude conductivity tensor determine what value of quan- 
tized Hall conductivity the system will acquire at large 
scales, when the interference effects suppress the diagonal 
conductivity. 

Basing on insensitivity to the character of disorder, 
we expect that weakly chiral network of junctions and 
point contacts represents electron gas in a non-quantizing 
magnetic field. To establish the correspondence we need 
to express the values of a^x and axy of our network in 
terms of parameters p and q in the Drude regime, i.e., the 
regime where the interference can be neglected. Above 
we have already identified {kpl)~^ with parameter p of 
the p-q model. Here we elaborate on this relation. 

Our main idea of modelling a competition be- 
tween interference-induced localization and magnetic- 
field-induced orbital curving is to confine the orbital ac- 
tion to a set of compact spatial regions-junctions, with 
asymmetry of scattering to the left and to the right pro- 
portional to magnetic field. In zero field, these junctions 
represent strong scatterers of electrons. These strong 
scatterers come in addition to weak "intrinsic" scatter- 
ers, that limit the mean free path, I, of electron gas. If 
strong scatterers are sparse, they will not affect the trans- 
port of electron gas at all. On the other hand, if they are 
dense, then the mobility will be limited exclusively by 
scattering off these strong scatterers. This suggests that 
the distance between the scatterers should be chosen of 
the order of I. Therefore, as we replace realistic electron 
motion by a motion along the links of a network, the di- 
mensionless lattice constant should be chosen as ~ kpl. 

In our p-q network Fig. [Sj in addition to junctions 
(strong scatterers) there are point contacts on the links. 
The role that these point contacts play is the following. 
Our junctions Fig. [S] do not provide any backscatter- 
ing. More precisely, as can be seen from Fig. [5l an 
electron starting along a given link, say, to the right, af- 
ter several scatterings off the junctions, will never return 
to the starting point from the left. Thus, the junctions 
alone cannot model the interference effects in realistic 
electron gas, where the probability of such a return is 
~ l/{kpl). It is the point contacts that provide possibil- 
ity of backscattering in the p-q network, and parameter 
p is chosen l/{kpl) in order to model the realistic return 
probability. 

From the above reasoning, the diagonal conductivity 
of the p-q model is axx ^ In a finite magnetic field, 
estimate for axy can be found from Eq. ([29]): 

(Jxy = RhCTxx ^ ( o ^ ?! '^xx ^ — -■ (72) 



Note now, that the boundary of delocalization transition 
established in the present paper, is a(l/2 — q) = p with 
a ^ 1. Then we conclude from Eq. (|72p that transition 
occurs when the "Boltzmann" value of axy is ~ 1. On the 
other hand, levitation scenario is based on the conjecture 
that the classical axy — 1 does not get renormalized upon 
increasing the sample size, see Eqs. ([S]), ((B]). Consistency 
between delocalization condition within p-q model and 
within scaling theory^ can be viewed as evidence that 
the p-q model adequately captures microscopic physics 
behind the levitation scenario. 

There is a fundamental reason why the picture of the 
weak-field quantum Hall transition is much more complex 
than the picture of the strong-field transition. Namely, 
for strong-field transition there is an exact duality of elec- 
tron states above and below critical energy. By con- 
trast, there is no such inherent duality in the weak- 
field transition. This, in particular, does not allow to 
employ the quantum real-space renormalization group 
approach2Si2&"-i2£ to describe this transition analytically. 

Duality with respect to the center, e = 0, of the Lan- 
dau level in the CC model insures that, with a strong 
spread in the saddle-point heights, the scaling region nar- 
rows, but delocalizcd state remains at e = 0. A remark- 
able outcome of our numerics is that the same property 
holds for the weak-field transition: we have verified that, 
upon introducing a strong spread in the local values of 
the backscattering strength, p, but keeping (p) fixed, does 
not effect the position, Pc{q), of the transition point. On 
the other hand, upon increasing the spread in the local 
values of p, interference effects become progressively less 
relevant. This allowed us to uncover a transparent clas- 
sical picture of the low-field quantum Hall transition, see 
Figs, m [ini which is a counterpart of the percolation 
picture of the strong-field transitioniSr^. The fact that 
levitation emerges both within the p-q model and from 
the scaling equations Eqs. dS]), ©, still does not mean 
that the p-q model offers a microscopic support for the 
Pruisken theory. After all, phenomenon of levitation fol- 
lows from a general argument that delocalized states have 
nowhere to go but up in the vanishing magnetic field. In 
fact, scaling equations Eqs. JS]), ^ suggest a stronger 
message: namely, positions of delocalized states do not 
change if magnetic field does not exercise any "phase" 
action. Neglecting the phase action corresponds to re- 
placement of the first term in Eq. ([5]) by a constant. We 
emphasize that the p-q model supports this prediction of 
the Pruisken theory. This is reflected in the coincidence 
(within accuracy of our numerics) of the p-q boundary 
with and without TRS. 

In closing, we list the questions which were previously 
addressed in the frame of fully chiral CC model and are 
also pertinent to the weakly chiral p-q model: 

(i). Two-channel network m odel s of CC type 

were studied in Refs. |2i|2l7l|l7 ll0ll - [T03l in connec- 
tion with spin-orbit induced splitting of quantum Hall 
transitional^, disorder-induced "attraction" of delocal- 
ized states from different Landau levels^S'^'iSi and quan- 
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turn Hall transition in disordered superconductor&i22iii2^. 
All these considerations are based on two co-propagating 
channels on each link. In p-q model, the two channels 
are counter-propagating and spinless. It would be in- 
teresting to incorporate spin-orbit coupling into the p-q 
model for the following reason. In zero magnetic field and 
in the presence of spin-orbit coupling there is a critical 
energy above which electron states are delocalized^iiSi. 
On the other hand, in strong magnetic fields, spin-orbit 
coupling splits discrete delocalized states^i^^. Upon de- 
creasing the magnetic field, both delocalized states are 
most likely to head towards zero-field metal-insulator 
transition point. Microscopic model describing this sce- 
nario must include orbital action of magnetic field, spin- 
orbit coupling, and interference effects. This can be ac- 
complished upon incorporating spin degree of freedom 
into the p-q model. A possible application of the above 
physics is graphene in a weak magnetic filediS^. In the 
latter case, the inter-valley scattering plays the role of 
spin-orbit coupling. 

[ii). In three dimensional layered system in zero mag- 
netic field, arbitrarily weak coupling between the layers 
delocalizes electron states above a certain critical energy. 
In strong magnetic field, a weak interlayer coupling smear 
discrete delocalized state into metallic band^"^. Match- 
ing of these two scenarios takes place in weak magnetic 
fields. The p-q model based description can be employed 
in this domain. 

{Hi). Interaction- induced dephasing^*'^ is crucial for 
experimental observability of levitation, since the transi- 
tion is smeared when localization radius exceeds the de- 
phasing length. The issue closely related to dephasing is 
a peculiar behavior of the non-diagonal resistivity, pxy , in 
the vicinity of the high-field transitioniSSrJJ^. The ques- 
tion about behavior of pxy near the low-field transition 
point can be asked and addressed within "incoherent" 
p-q model. In addition, the question about interplay of 
disorder and interaction previously studied for high-field 
transitionii^"-i2^, can be redirected to the low-field tran- 
sition. 

{iv). It is known that there is an intimate relation be- 
tween delocalization transition in CC model and critical 
behavior of superspin chainsi2^"-i2^. It would be interest- 
ing to investigate wether the p-q model corresponds to 
any spin model. It is also known that CC model is re- 
lated to the Dirac Hamiltonian with disorder-i22ii^. Re- 
call that at degeneracy point, p = 0, g = 1/2, the p-q 
model falls into two independent CC models. Modifica- 
tion of the mixing of two counter-propagating channels 
on the links transforms the p-q model with q = 1/2 to 
the model studied in Ref. 59. This modification leads 
to criticality which was studied within the sigma-model 
approachl^iii^. It would be interesting to investigate 
such a modification of the p-q model away from the de- 
generacy point. More specifically, the case q ^ 1/2, in 
our p-q model, differs from the model Ref. |59| primarily 
in the scattering matrix of junctions. Our junction ma- 
trix Eq. (1^51) describes scattering, say, to the right, with 



the same probab ility q for all incoming channels, while 
the nodes in Ref. [S^ describe scattering to the right with 
the probability q from two opposite channels, and with 
the probability (1 — q) from the other two opposite chan- 
nels. In terms of coupling of counter-propagating chan- 
nels on the links, in Ref. [59| this coupling differs from 
our link matrix in following respects: to get the cou- 
pling of Ref. 59 one should replace in Eq. (1551) the phases 
should all be put zero, while satisfy the relation 
<y5i — '^2, ¥^3 — which is different even from the case 
with TRS in the p-q model. 
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X. APPENDIX 

The easiest way to derive Eq. (IT^ is to notice that 
with respect to the motion along the diagonal direction of 
network the diffusion has, effectively, a one-dimensional 
character-i2^. This is because the motion along the diago- 
nal can be viewed as a random sequence of transmissions 
and reflections without deflections. For a given site, the 
full probability of transmission along the diagonal direc- 
tion (say, up and to the right) \s T = t^ -\- (P , while the 
full probability of reflection is 1 — T = + . Parameter 
kpl should be identified with ratio, T/(l — T). Together 
with flux-conservation condition, + + 2(P = 1, this 
ratio reduces to Eq. 

As we pointed out in Section II, the point d = l/-\/2 is 
singular: strong localization predicted by scaling theory 
with Boltzmann result Eq. (jl2p as an initial condition, 
contradicts to the result of Ref. i59|. This contradiction 
can be resolved from the following reasoning. 

The same average {(P) can be realized in two com- 
pletely different ways: 

(i) Boltzmann approach. In the two instances of scat- 
tering at the same node, electron can be deflected both 
to the left and to the right with equal probabilities, 50%. 

(ii) Percolative approach. A given node scatters always 
to the left or always to the right] 50% of the nodes scatter 
to the left, and remaining 50% scatter to the right. 

Note, that the classical motion of electron is very dif- 
ferent in the above two realizations. In the case (?) this 
motion is diffusive. Then interference will readily localize 
this motion. The case (iz) corresponds to classical per- 
colation threshold. To see this, notice, that the electron 
predominantly travels in loops (clusters) and one trajec- 
tory is the inflnite cluster spanning the entire system. In 
fact, it can be easily seen that possible trajectories belong 
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to two decoupled bond-percolation networks. Interfer- 
ence will transform these networks to two CC models at 
the delocalization-transition points. Since the percolative 
picture is fully coherent, it captures properly the quan- 
tum delocalization, while "dephased" Boltzmann picture 
does not. 

Thus the distinguished characteristics of the point d = 
I/V2 is that only at this point there is sensitivity to the 
realization of the disorder. 



Note that in the case t = r = but di = y^l — d"^ not 
equal to l/-\/2, prediction of the scaling theory that the 
Boltzmann result governs localization properties at large 
distances is correct. Although with a strong spread in 
local d- values, when the network still breaks into two CC 
networks, the states of both CC networks are localized 
with localization radius, ^ \df — l/2\~'^/^. 
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